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ABSTRACT 


A computer-controlled system for automated, remote 
measurement of phase and amplitude shaded doublet, triplet, 
and quadruplet acoustic arrays beam patterns was developed 
and tested. Experimental data were collected with an 
anechoic chamber to minimize backscatter and other 
interference effects. 

Preliminary impedance, reactance and resistance plots 
were made on several speakers to ensure that a matched set 
with similar characteristics was utilized for the arrays. 
Once the arrays were established, a computer program was 
developed for an HP-86 computer which allowed it to rotate 
the arrayS, sample and collect data, and plot the results. 

Various combinations of amplitude and phase shading were 
applied to the dipole, tripole, linear quadrupole and 
quadratic quadrupole array. The results were then compared 
to computer generated beam patterns which agreed to + 2dB. 
This uncertainty was mainly generated by the slight 


deviations in the positioning of elements. 
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A. PURPOSE OF STUDY 

The purpose of this investigation was to develop and 
verify a computer-controlled experimental system to measure 
and plot the pressure fields of various acoustic arrays in 
an anechoic environment. This was to be done using various 
combinations of phase shading and amplitude shading of 
doublets, triplets, and quadruplets. Experimentation was 
conducted under computer control to maximize the amount of 
data available and to ensure the highest degree of accuracy. 
Additionally, a comparison was conducted with the 


theoretical predictions obtained by Park [Ref. 1]. 


B. BACKGROUND 

The use of dipoles and quadrupoles for beam forming and 
direction finding, both in air and in water, has been a 
Subject of intense investigation since the 1948's. The 
basic radio direction finding apparatus of the 1948's was 
nothing more than a set of orthogonal dipoles’ which 
consisted of two crossed loops. The relative signal 
amplitude from the two dipoles was used to determine the 
direction of arrival of a signal. The angle is equal to 
Arctan (y/x) where x was the amplitude of one dipole and jy, 


e 
the amplitude of the other dipole. This presented a 188 


im 


bearing ambiguity which was resolved by comparing the phase 
of the x and y dipole signals with the phase of a co-located 
omni-directional antenna. [Ref. 2] 

In the 1958's, extensive studies were conducted using 
planar arrays with multiple acoustic and electromagnetic 
dipole elements. These studies led to the formulation of 
planar arrays employed in radar systems of the 1969's. 
Further research resulted in the ability to steer beam 
patterns by use of multiple elements, phase shading, and 
amplitude shading. 

In 1976, work was done by Hughes and Thompson to 
investigate acoustic communications, particularly in 
underwater acoustics for cases where the position of the 
source or the receiver is fixed. They investigated the 
tilted directional response patterns below the horizontal 
plane due to amplitude weighting and a_ single 99° phase 
shift. Here, a method of achieving a tilted pattern from a 
symmetrically shaped array without using numerous’ phase 
shifts or time delays was discussed. [Ref. 3] Moses and 
Smith also examined beamforming this time with a cylindrical 
array uSing real (rather than complex) shading coefficients 
[Ref. 4]. 

In 1985, Flanagan examined the use of computer-steered 
microphone arrays which were two-dimensional planar arrays 
of 9 x 7 elements. Here an array was designed and computer 
steered to discriminate against sound arrivals from all but 


the desired directions. [Ref. 5] 


One example of the Navy's current use of dipole 
technology is in sonobuoys. Hive SSQO-53 DIFAR Sonobuoy 
possesses both an omnidirectional and bi-directional 
hydrophone capability. The SSQ-77A VLAD Sonobuoy possesses 
an upper and lower array strung ina single line array. 
Beamforming is accomplished by the use of both phase shading 
and amplitude shading. 

These advancements in the study of arrays ranging from a 
Single dipole to large scale planar arrays, coupled with 
computer control, form the basis for this further study of 
shaded arrays. There is currently a related research effort 
Ongoing in the study of computer controlled arrays and beam 
patterns being conducted by the Environmental PhySics group 


at the Naval Postgraduate School. 


C. REPORT OVERVIEW 

In Chapter II, a basic introduction to theoretical 
pressure and directivity calculations of a dipole is 
developed. This is then used as a base to expand to the 
physical case and the electrical linear systems case. 
Finally, these results are applied to different forms of 
quadruplets and triplets to show theoretically, the 
expected results. 

Chapter III contains a summary of the beam pattern 
measuring system, the computer program development, 


experimental procedure, and representative data and plots. 
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Additionally, it provides a description of the system which 
was assembled for measuring the pressure fields of the 
various arrays. 

Chapter IV contains a summary of data analysis and 
verification. Additionally, representative plots are 
provided for amplification of arguments. 

Chapter V contains concluding remarks and 
recommendations for further work. 

Appendix A provides a brief description of each piece of 
equipment and summary of its capabilities as applied to this 
experimental investigation. 

Appendix B lists the BASIC computer language program 
used for data acquisition and formulation in this study. 

Appendix C contains specific plots, both experimental 
and theoretical, obtained as a result of this investigation 
and that of Park [Ref. ll]. These plots were obtained for 


beam patterns of various array configurations. 
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LILA eeWclO log 


The dipole has been extensively explored and is now 
utilized as a basic unit in the construction of a great many 
arrays and antennae. In this section we will elucidate the 
theory using two approaches 1) an acoustic approach in 2 
dimensions and 2) a linear systems approach in three 
dimensions. These two approaches are related to each other 
and yield the same result after minor changes are made in 
the geometry of the problem. Both are presented because of 
direct Navy applications dealing deal with SONAR arrays. 

Following the presentation of each of these, the 
equations for the triplet and quadruplet arrays utilized in 


this research effort will be shown. 


A. DIPOLE RADIATION-ACOUSTIC APPROACH 

Because it iS usually desirable and useful to be able to 
determine the direction from which an acoustic’ signal 
Originated, arrays were developed to provide this 
information; a simple example of this is the acoustic 
dipole. A simple dipole consists of two simple sources of 
equal strength but opposite phase, separated by a distance d 


pelt (2. 1)". 
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Figure 2.1 Geometry used in Deriving the 
Radiation Characteristics of an 
Acoustic Dipole. 


The pressure, p(r,@,t), measured at the field point (r,@) 
due to source 1 and 2 respectively are 

Pr = {A(r + Ar, )] expfi(at - k(r + Ar,)] (1) 

Py = [-A‘(r - Ar,)] expfi(a@t - k(r + Ary) (2) 
where A is the pressure amplitude of an individual source 
measured at one meter from the source, r is the distance 
from the field point to the midpoint between the sources, 
W=2 Jif is the angular frequency and k is the wave number. 


O 
The minus sign associated with P2 accounts for the 1898 


phase difference. The total acoustic pressure is the sum of 


L6 


Evmeana P2, 


P(r,8,t)=(A‘r){{fexp(-jkAr, )] [1+(Ar, /r)]} - {(fexp(-jkAr,)] [19 (Ar,/r)]} }exp(j(@t- kr)) (3) 


Considering only field measurements where the separation d is 


small compared to the distance r, we can use the far field 


, ) 4 y | 
approximations, Ar, , = (d/2)sin8 and An git << : 


This yields the acoustic dipole radiation pressure field, 


P(r,8.t) = -j (2A/r) sin {(kd.2)sin 8)} exp j(@t - kr) (4) 


- 180 -30_— 0 90 180 #8 
(a) 


Figure 2.2 Cartesian and Polar Representation 
of Pressure Amplitude Distribution 
of Dipole Radiation. [Ref. 6] 
VW 


As noted in Eq. (4) and Seen in Fig. (2:2) )plottiunaqeenemnenn 
field for d<<A, pressure nulls exist when (kd/2)sinQ =7J/n 
where n is an integer. [Ref. 6] 

The twin lobes and nulls of a dipole result in an 18 
ambiguity when trying to localize a source. If the dipole 
is combined with an omnidirectional receiver, the ambiguity 
can be resolved. In Fig. (2.3), various combinations wer 
dipole and omnidirectional receivers are combined to provide 


a cardioid response. 


BIDIRECTIONAL 


NONCIREC TIONAL 





COMBINATION 





Figure 2.3 Directional Diagrams of Various Combinations 
of Bidirectional and Nondirectional 
Microphones. [Ref. 7] 


Heaney is the voltage output of the omnidirectional 
receiver and Ry the voltage output Of the directional diver 
recelver. It is not important for the two sources to be of 
the same sensitivity, but the ratio R,/R2 must be unity for 


true cardioid response. The phase shift between these _ two 
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receivers, which yields directionality, is achieved through 
a difference in the phase characteristics of the omni which 
is differentiated (high pass filtered) to compensate for the 
linear frequency dependence of the dipole sensitivity. [Ref. 
7). 

This design is presently used in SSQ-77A VLAD Sonobuoys 
where a cardioid is formed pointing in the downward 
direction. This maintains the null in the upward direction 
in order to alleviate high frequency interference from wind 


and sea surface noise. 


Be LINEAR SYSTEMS 
Consider the geometry for a dipole as shown in Fig. 


(2254): 





Figure 2.4 Dipole Geometry [Ref. 8] 


arg 


In this representation as developed in Ref. (7), C) andes, 
are complex weighted point sources which are free to rotate 
in the x-y plane and d is the distance between elements. 


The complex weighted response factor ci Of the ith point 


source is 
Ci n a. exp(j@. ) (5) 
where ai is the amplitude and Wea is the phase angle. For a 


dipole, the complex weighted response factors are Cie 
and ‘%.; = -a, « 

For a wave originating from a distant source and 
propagating to the receiver, it is observed that there is an 
essentially constant pressure magnitude over the dipole 
aperture. Because the phase varies with position along the 
wave front in a traveling wave, the response of the dipole 
is given as a function of frequency and direction cosine u 
where u = sin@cosY. The response is 

Diu) 1a, exp(jtdu/A) + c, exp(-jmdu, 4) (6) 

This equation is simplified for the dipole to: 
D(ffu) = 2ja, sin (mdu; A) (7) 
Since u is a direction cosine, it is only valid over the 
interval ~l<u<l. As a result, 0SO0Snx and OSWS2% , y=o 
implies a broadside or on-axis geometry and u = 1 implies 
an end-fire geometry. A plot of D(f,u) as a function of u 
can be used to approximate the response of a dipole. [Ref. 


8 ] 
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An example of this two element response of two point 
sources with a separation d equal to one half the wavelength 
of the traveling wave results in a real response for values 
of u from (-A/2d)to(+A/2d) - If the sources are assigned 
dipole weightings, substitution of the functions © and Y for 
direction cosine u result in the unnormalized, far field 


beam pattern 


D(fO.ay) = j2a, sin[(md/d) sin 8 cos w] (8) 


This linear systems theory can be easily applied to the 
axial quadrupole and Tesseral quadrupole. The axial 
quadrupole is an array of three complex weighted point 
sources where Cc =a,=c,=a,=1 Aa) Ge = ay = =) ae ip 
fg, (2.5). 


Here Co can also be thought of as 2 sources with each having 


an amplitude weight equal to -l. 


(2,9, w) 


Figure 2.5 Axial Quadrupole Geometry [Ref. 8] 


Jad 


With the 1,-2,1 amplitude weighting, the far -£ lelar 
unnormalized directivity function is defined as: 
D(Ef.) = [exp(-j2f,d'2)] - 2 + [expj2nt,d: 2)] (9) 
where f, = uA = sin@cosy A , This reduces to: 
D(f6.w) = 2 [cosi(md A) sinBcosy} - 1] (18) 
In the case of a Tesseral quadrupole Fig. (2.6), a 
planar array of four amplitude-weighted point sources, the 
unnormalized directivity function is: 
D(ff fy) = 2Q cos[md(f, “ a) - 2Q cos[md(f : 9) (11) 
where f= wA= sin@coswA . 6 f = wd = sinOsinyA and the 
direction cosine v = sin@ sinY. This yields 


Dif O.w) = -4Q sin[(md'A) sinBcosy] sin{(md/A) sin@sinwy] (ik) 


(2,9, ) 





Figure 2.6 Tesseral Quadrupole Geometry 


a 


which is the far field, unnormalized directivity function 
for this quadrupole as a function of @ and Y, (Ref. 8] The 
above development 1S particularly important because this 


type array was investigated in the course of this work. 


C. GENERAL QUADRUPOLE 

Park, in his thesis (Ref. 1), developed an equation for 
calculation of a four element array Fig. (2.7). The 
equation was developed by adding complex pressures from four 


sources to yield 
p = ze el(wt-kKR-kD,sinOsin®) , me el(wt- KR + kD,sinB@sin®-@,) (73) 
R R : 


h: - oi(wt -kR + kD,sinOcos® - @,) 


f ' el(wt - KR - kD sin@cos® - @,) 


where [(1+A+B+C) R] exp{j(mt - KR) is the maximum possible 
pressure amplitude and the distances Dj, Do, D3, Dg and 
measured from the origins to the particular element. The 
Gueeectivity is: 


I ae ane (14) 
H(0,D) = ——__—__——_{[cos(kD, sinOsin®) + Acos(kD,sin@sin® - @,) 
een 2 2 


+ Bcos(kD ,sinOcos® -@;) + Ccos(kD ,sinOcos + QP 
+ [sin(kD sinOsin® ) . Asin(kKD,sinOsin'P - @,) 
- Bsin(kD ,sin@cos - @,) + Csin(kD sinOcosD + @ 7} ‘. 
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Figure 2.7 Quadrupole Geometry 


If Dy, = D9 = D3 = Dg then comparison of Figs. (2.7) and 
(2.6) shows that the configuration is simply the Tesseral 
quadrupole of Fig. (2.6) where V has been shifted Asan 
Additionally, the axial quadrupole of Fig. (2.5) is 
obtained using Park's program by placing source B at an 
infinitesimally small distance from the array center 
(i.e., D3 = 6.00081 meters for any A) with amplitude 
weighting equal to twice sources 1 and A and _ source C 


amplitude weighted to zero, which removes its contribution 
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from the array. This axial quadrupole can also be obtained 
by placing sources B and C at the array center and amplitude 
weighting equally. 

These equations were utilized to provide verification 
data for the laboratory experiments. The three element 
array data were verified by amplitude weighting element C of 
Eq. (14) to zero. Dipole data were verified by amplitude 


weighting elements B and C of Eq. (14) to zero. 
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III. EXPERIMENTAL PROCEDURE AND DEVELOPMENT 


A. BEAM PATTERN MEASURING SYSTEM OVERVIEW 

The beam pattern meaSuring system is an instrumentation 
package which was designed to power and rotate an acoustic 
array and then receive the signal, compute the pressure 
distribution of the array and display this data in both a 
tabulated format and plot of sound pressure level in 
decibels versus degrees of rotation. All data nege taken in 
an anechoic chamber to minimize surface interference. 

The system used for meaSuring beam patterns of different 
arrayS is divided into two major groups. These are the 
Signal generation and input section and the automated 
advancement, data collection and display section. An 
overview of these is now presented along with a block 
diagram of the system, Fig. (3.1). A detailed description 
of each piece of equipment is contained in Appendix A. 

The Hewlett-Packard (HP) Function Generator Model 3314A 
was utilized for all signal generation. The signal passed 
through an individual HP-467A Power Amplifier to each of the 
University Sound Model ID-38 acoustic driver units. Setting 
the output levels of the driver unitS wasS accomplished 
utilizing the variable gain control on the power amplifiers. 


This constituted the signal generation and input section. 
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Figure 3.1 System Diagram 
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The receiving and processing section was controlled by 
the HP-86B computer. It directly controlled the array 
rotation and the reading of all measurements of peripheral 
instruments via the Hewlett-Packard Interface Bus (HPIB). 
The computer sent clockwise stepping commands to the motor- 
power/direction controller via the HP-3421A Data 
Acquisition/Control (DAC) to rotate the array in stepped 
increments of about one degree. Following each advancement 
of the array, the computer read both the resistance output 
of the array rotation measurement potentiometer and the 
voltage output of the General Radio Model 1962 Microphone. 
The array rotation potentiometer was a @ to 18 Kohm 
potentiometer mounted on top of the rotation motor and 
measured the angle of rotation of the array (@ to 18 kKohms 
corresponded to 8° to 360°). These resistance and voltage 
readings were then stored in the computer for use in the 
data analysis and plotting routines. Once the data was 
taken for 360° of rotation, it was used to compute the 
rotation angle and its corresponding decibel (dB) level for 
each datum. These values were then output to the computer 
monitor and plotted in the form of a polar plot of dB versus 


degrees of rotation. 


B. BEAM PATTERN PROGRAM 
The Beam Pattern program listed in Appendix B, consists 


of a series of sub routines. The program starts by 
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providing addresses for the two multimeters and the data 
acquisition/control unit. Next the initialization 
Subroutine resets and initializes these three instruments. 

Once the equipment is ready to operate, the incrementing 
subroutine is called. This is an iterative and interactive 
routine used to rotate the array and measure the voltage of 
the sensed pressure field and the corresponding resistance 
reading which is proportional to the angular position of the 
array. Up to 368 data points can be collected per data run. 
After the initial datum is read, the array rotator motor is 
energized for 9.5 milliseconds to increment the array 
azimuth. The array is then allowed to coast down for 4 
seconds prior to reading the voltage and resistance. This 
selection of rotate and coast down times is discussed later 
in the chapter. Once the resistance and voltage readings 
are taken, they are stored in an array for later transfer to 
a data disk and for use in the plot subroutine. 

Finally the plot subroutine calculates the relative 
decibel level of the beam pattern based on the measured 
voltages. It then takes the resistance reading and converts 
it to the angle of rotation. The calculated decibel levels 
are then normalized to zero decibels based on the highest 
received voltage reading. Finally, the data points are 
plotted on a polar plot of relative sound pressure level in 


decibel versus angle of rotation. 


29 


All raw data is stored on a data disk after a complete 
data run. Each data run is stored by run number, as shown 
in. Table gig It can be retrieved and printed at the 
printer as received voltage and resistance for each datum 
point or printed and plotted as decibel versus angle of 
FOtatdon. 

The physical weight of an array was a major variable in 
this program with regard to the selection of the number of 
data points needed. Arrays weighing less than 4. pounds 
required 368 data points for a full 368° of rotation. The 
number of data points required for larger arrays was 
decreased as array weight increased. For the four element 
array, weight was substantial enough (19 pounds) to require 
only 317 data points for a full 360° rotation. This was 
caused by the fact that the heavier arrays would coast-down 
for a larger distance after the 9.5 millisecond motor-on 
period. 

Another contributor to variation in the required number 
of data points was array balance. When utilizing a three 
element array, the array was sufficiently out of balance to 
require a greater number of data points for full rotation. 
This was attributed to minor binding of the rotator shaft as 
the array coasted down. 

Selection of motor run time and coast-down time was 
based on experimental investigation. A 998.5 millisecond 


P 0 : 
motor on time provided 1 rotation regardless of coast-down 
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time with a light weight array. The 4 second coast-down was 
selected because smaller periods did not allow the array to 
stop and the resistance multimeter to stabilize prior to it 
being read. Motion at the variable potentiometer when the 
resistance values were read yielded erratic reading on the 
multimeter due to electromagnetic interference generated by 


the motor. 


C. EXPERIMENTAL APPARATUS DEVELOPMENT 

The experimental apparatus was designed to optimize a 
thorough search of two, three and four element arrays and to 
establish their respective beam patterns. Previously, 
several experimental investigations have been conducted 
utilizing a two microphone array with a single source. 1 = 
was, therefore, decided that reciprocity would be invoked 
and the array would consist of multiple ID-38 driver units 
(Appendix A) with a single microphone utilized for sensing 
the pressure field. 

Six ID-38 driver units were checked individually in the 
anechoic chamber for their beam pattern impedence, 
resistance and reactance characteristics. This data was 
utilized to select the best matched set of four driver 
units. All drivers were examined between 5 Hz and 30098 Hz 
in an effort to locate and identify the resonances in the 
individual speakers. Fig. (3.2) shows a _ representative 


resistance curve for these drivers between the above 
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frequencies. It was desired to conduct all analysis at 
frequencies in relatively flat areas away from resonances. 
This was done to assure phase stability since reactance is 
net Varying. The speakers utilized in this investigation 
had resonances as per Table I. 

TABLE I 


ID-38 DRIVER UNIT RESONANT FREQUENCIES 


peaker Number : Fj Resonance (H : Second Resonance (Hz 
ik 689 2290 
Z 669 2225 
3 788 2159 
4 798 ZNeZ.9 
=) 746 PPS 
6 778 2125 


Matched set number 1 consisted of speakers 1 and 2 and set 
number 2 of speakers 3 and 4. All array beam patterns were 
run at 388 Hz and 968 Hz which fell in the flat areas of the 
impedance curves. 

The array existed in two individual configurations, a 
linear and a crossed axis system with variable tilt (Fig. 
oe 3) ) The linear system was used for investigation of 
doublets, some triplets and quadruplets (in the form of an 
axial quadrupole). The crossed axis array proved to be 
extremely versatile and useful, providing four arms, each 


78 centimeters in length, for holding the drivers. 


og 


Additionally, it provided a means to vary the angle®@ a full 
99°, Each driver was equipped with a harness assembly which 
mounted to the array rod in such a manner that the driver 
face was located 2.5 centimeters off the axis. This did not 
prove to be a matter of concern since all pressure 


measurements were made in the far field. 





e 
y 
LA 
+ 





Quadruplet 


Figure 3.3 Array Configurations 
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A scattering and shading problem of small magnitude was 
encountered because of the large center piece on the crossed 
axis array. This scattering problem was largely resolved by 
wrapping the center piece with fiberglass. 

The General Radio Model 1962 microphone was located 3.2 
meters from the center of the array placing it in the far 
field. It wasS mounted at the same height in the anechoic 
chamber as the array center. Changes in tilt angle theta 
were accomplished by tilting the array. 

The outputs of the ID-39 driver units were balanced 
prior to each data run. It was done by aligning the array 
in the vertical plane with the drivers facing the 
microphone. Then the output of each driver was adjusted 
utilizing the variable gain control on the power amplifier. 
This allowed all the drivers to be driven within a received 
voltage of 9.982 VAC of each other out of a typical level of 
i Volt. 

Output signals from the HP-467A Power Amplifiers were 
monitored on an oscilloscope as was the received signal from 
the microphone. Because the Signals were clean and without 
any substantial noise, filters were not utilized in the 


microphone circuit. 


IV. DATA ANALYSIS A VERIFICAT 


After establishing the ability to accurately rotate an 
acoustic array while collecting and displaying sensed 
pressure levels under computer control, it was necessary to 
ensure the data were accurate. This could be done by taking 
a set of parameters and obtaining both theoretical and 
experimental beam patterns for them. Verification was 
conducted utilizing the results of a FORTRAN program 
developed by Park (Ref. 1) which yields a variety of useful 


plots with the help of DISSPLA. 


A. THEORETICAL BEAM PATTERN PREDICTION 

With the program developed by Park, quick beam pattern 
plots could be made to compare with experimental patterns. 
As an example, Fig. (4.1) depicts a three dimensional 
representation of a dipole. Clearly, this particular is not 
of much use. 

The three dimensional plot was sectored for greater 
ues lity. Fig. (4.2) shows a sector representation for a 
four element array with kD) = kD) = 3.52 and kD3 = kD, = 0. 
Note that the kD; utilized in the remainder of the thesis 
related directly to the spacing from the center of the array 
to the element as related to Fig. (2.7). Refer to Fig. 


(2a The Z-axis of the plot is beam pattern (dB), the Y- 
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Figure 4.1 Three Dimensional Dipole Representation 


og 



































ANIA gz 

; ay iia HEE: 
Xe | 

a La We ri f PR 
a My fi Aas Uh 
ES sy \ AF * 
Soy Han W 4 
aS Bs 4 | , 





W A 


ny a 
‘i Wh 


ey 
i 


kD1=kKD5=kD3=kD,= 10.55 

DI=06 D2-06 

D3= 0.6 D4= O. 6 

VIEW ANGLE=+20.+60 

er O P3=0 P4=9 
A=1. B=1. C=1. 


C 
Figure 4.2 Sectored Three Dimensional Plot ( 4 Element ) 


38 


axis psi (degrees), and the X-axis’ theta (degrees). 
Directivity is scaled from 28 to -48 dB (with 9 dB being the 
highest achievable value because of normalization), psi from 
180° to 360°, and theta from 0° to -99°. This plot appears 
to be very busy but proves to be quite useful. It is easy 
to evaluate this sectioned representation by following the 
contour line of the directivity above the line theta equals 
9° for psi from 180° to 369°. This provides the same dB 
levels at 189° < psi < 360° as the polar plot representation 


plotted in Fig. (4.3). 


B. EXPERIMENTAL DATA-SINGLE ELEMENT DIRECTIONALITY 

As stated in Chapter III, testing was conducted at 389 
Hz and 969 Hz which were in the flat areas of the drivers' 
resistance curves, away from their resonance peaks. This 
was done to assure phase stability between drivers since the 
reactance 1s not varying (large phase shifts do not occur 
except near resonance). 

Initial data runs were conducted with a single ID-39 
acoustic driver to establish the directionality of each of 
the drivers. The drivers were rotated 360° at three 
separate tilt angles of the driver face, @ = a°, 45° and 90° 
A representative graph of the directionality is contained 
in Fig. (4.4) which was made with the driver at a tilt angle 
of 9° and pointing at the microphone prior to rotation. The 


dotted line is the laboratory beam pattern. Note that there 
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is a l to 2 dB drop below the O dB level as the driver 
points away from the receiver which shows the driver is not 
completely omnidirectional. It should also be noted that 
the program's normalization routine sets’ the highest 
received microphone voltage level to 9 dB and normalizes all 
other voltage levels in the beam pattern to this level. 
This was done to provide greater ease in plotting of beam 
patterns so as to produce a uniformity between plots. 

C. EXPERIMENTAL AND THEORETICAL DATA - VARIOUS ARRAY 

CONFIGURATIONS 

An investigation of two, three and four element arrays 
was conducted. For each of these configurations, several 
different combinations of spacing and amplitude and phase 
shading were looked at. Table 2 at the end of this chapter, 
is alist of these variations for which beam patterns were 
obtained. Note that phase shading was limited to either g 
or 189°. In viewing the figures, the term "x-pole" is used 
loosely to represent the term “x-let" (i.e., quadrupole for 
guadruplet). 

For the two element array, beam patterns were run in the 
laboratory for both- in-phase and 188° out-of-phase doublets. 
An example of the out-of-phase doublet beam pattern obtained 
in the anechoic chamber is shown in Fiq. (4.5). (Note that 
in this figure kd is the quantity 2kD, when Dj = Do.) 

Its theoretical verification, obtained from Park's program 


is shown in Fig. (4.6). A contributing factor to their 
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small differences, is the slight directionality associated 
with the acoustic driver. The slight offset of the 
laboratory plot nulls is caused by the array not being 
exactly aligned to broadside at 8 Kohm on the potentiometer. 
The 99° rotation between the two plots is caused when theta 
is fixed at 98° vice 9° for theoretical plots. Because of 
the narrowness of the nulls, laboratory data would require a 
sample exactly with the array broadside to achieve the null. 
Other than the difference in nulls, there is exceptionally 
good agreement between experimental and theoretical beam 
patterns. 

In the investigation of three element arrays, the arrays 
were configured as equilateral triangles, isosceles 
triangles and as a linear array. For phase shading, three 
cases were looked at: all elements in phase, the center 
element 189 out- of-phase, and one end element placed 189 
out-of-phase. These arrays were also subjected to various 
amplitude shading (the axial quadruplet composed of three 
elements with the center element amplitude shaded to twice 
the value of each of the end elements will be discussed in 
the quadruplet section). Experimental data for a three 
element isosceles triangle with all elements in phase and 
KD) = kDy = 2kD3 = 5.28 are shown in Fig. (4.7) with its 
theoretical verification shown in Fig. (4.8). The deep 
theoretical nulls seen on Fig. (4.8) are not present in the 


laboratory data. This can be attributed to a small variance 


45 


JOTd SOTSOSOST JuaWeTY ser, AroQeIOGeT /*p eanbtTy 


Beg g=379NY 2°/A SII ‘E =TIdWy 5289 0 'O ‘O=3SVHd 
Giz°s ="quz="aqyealqy! WwW O2°E=4 





CHP) NdsllVd WVsad 3 lQqalal 


46 


THETA VS. DIRECTIVITY 
0 





in the equal amplitude shading, a small amount of scatter 
from the large array center joint, and a Slight 
directionality of the individual acoustic drivers. 

The beam patterns of atriplet in the form of an 
isosceles triangle are shown in Figs. (4.9), (4.18), and 
(4,11) to show the effects of tilting the array out of the 
horizontal to angles of theta equal to 30° ; 60° and of 
respectively. The plane of the elements was initially 
perpendicular to the x-y plane, making theta, the angle 
between the z-axis and the array plane equal to 9° The 
array plane was then tilted 39° about the x-axis and rotated 
360° about the z-axis to obtain the data for Fig. (4.9). 


The array plane was then tilted in 38 increments to obtain 


Figs. (4.18) and (4.11). For this triplet case, kD, = kD, = 


2kD3 = 5.28 with kD3 sensing the elevation change (Fig. 
(45 ae) 


48 


~ 


eJOUL UT ITTL (O€ FOF UASIIed ureog 6°? sanbry 


Beq C€=3T9NY 9PA 20¥°2 =TdnV 5eq O‘O ‘O=3SVHd 
9iz7°s =tayz=Cay=lay = COZ "EU 


(dP) NdallVd 





WY¥3dG 4 lOdldl 


49 


PFOUL UT FTTL (09 OF uUTSIIed weSg OTP Sanbt4 


Bog gg“379N¥Y 2°A 298°! =TdHv Seq o'0 ‘D=aSVHd 
o1z7°s =“ayz=Cay=Tay sw COO2 *E=8 


——— — 


— 
“ 


a © aed — . a 
es vA > - se . \ : \ 
i: y 7 | 7 ~ . | 
/ i ye a ‘ \ \ 
on ae \ 
3xTH [ff ~ — } >: ema ee ae ) | 
ae. GOs 
\ \ . io ial fy | 
\ \ ; / H 
\ \ s 7 
\ \ ~ ea ae vs if 
. \ | ieee 
v “. ~ To 
\ = ; ee 
\ aa -s 4 
S | “ 
= a “ “ 


Cay) Neeeeeye VWeaie = etal 


20 


eCVDUL UT ATUL 906 AOJ urS}Weg Weeg [I’p eamnbty 


68q O6=379NV 2°A GB2°1 =IdWv 52g 0'0 ‘0=3SVHd 
giz7°s¢ =tanz=ca= lay H OC °E=¥ 





(dP) NdsllVd WV3a 3 lOalal 


Bal 





Figure 4.12 Array Configuration for Tilt Theta 


These variations in tilt angle (theta) show that the 3g 
tilt has a diminished dB level for all lobes except the 
broad lobe coming into and exiting broadside. This is due, 
in large part, to the directionality of the drivers, because 
they are facing the microphone to a greater extent at 30 
than at 99° of tilt and mutual shading exists. 
Additionally, for these various tilt angles, the Z-axis is 
rotating about the local vertical, as explained in section D 
of this chapter. 

The last major configuration investigated was the 
quadruplet. It was looked at both as an axial quadruplet, 
Fig. (2.5), and as a planar array, Fig. (2.7). An 


experimental beam pattern for a planar array quadruplet with 


KD) = kD2 = kD3 = kD4 = 5.26 and all elements in phase is 


a2 


Shown in Fig. (4.13) along with a theoretical beam pattern, 
Fig. (4.14). There is good agreement between the two plots. 

A plot of an axial quadruplet represented by an 
amplitude shaded triplet, obtained in the laboratory is 
P@ewnin Fig. (4.15). Comparing Fig. (4.15) to Fig. (4.13), 
it is seen that when two elements of the planar array are in 
line, the arrival from the more distant source is weakened 
by the shadowing effect of the nearer source. Therefore, the 
near source can be emphasised by as much as one dB-~ when 
compared to the far source. Additional plots of various 
configurations are contained in Appendix C, where _ good 


agreement continues to be seen. 


D. EFFECTS OF ARRAY TILT, THETA 

Referring back to Fig. (4.2), it is again noted that 
contour lines are directly related to polar plots’ when 
either theta or psi is fixed. However, the countour 
representation is not as useful when theta and psi are both 
changing with rotation of the array as in the case where the 
array has been tilted in theta. In this case the polar 
curve is found by superimposing a trajectory which is’ the 
specific locus of @ and Yin the © -\Y plane onto the 
contoured surface. 

For this work, the microphone was fixed in the far field 
at a distance of 3.2 meters and level with the center joint 


of the array. Tilting the array in angle theta caused an 
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offset between the Z-axis and the local vertical about which 
the array rotated. This can be seen in Fig. (4.16). In the 
left sketch, the array elements are aligned in the standard 
rectangular coordinate system, with the elements in the x-y 
plane and rotating about the z-axis. However, when the 
plane of the array is tilted 30° about the Ke ax Suey = v+30 
and z' = z+3@. A 96° rotation about the z-axis to 278 
yields z' = 98° - 30° = 60° . This means that a 360 
rotation about the z-axis causes Yl and Zl to vary from 126 
to 68° and back to 120° . ‘This is due to the fact that the 


microphone is fixed in space in the far-field rather’ than 


tilting with the plane of the array. 


a7 





Figure 4.16 Z-axis Offset Effects 


This problem does not occur in the classical representation 
where the microphone is fixed with respect to the axis of 


rotation of the array. 


Boo “BRRGRS 

The sources of errors in the beam pattern measuring 
system were examined both mathematically and empirically. 
Sources of error for direct error analysis were based on 
multimeter values and power amplifier settings. With output 
reading from the multimeters and power amplifier exhibiting 
precision to the third decimal place, instrumentation errors 
from these sources are less than 1% and are not major 
contributing errors to measured beam patterns. 


58 


The major source of error in the system was determined 
empirically to be the acoustic driver positioning on the 
array. The effects of a small change in kd are shown in 
Bags. (4.17), (4.18), and (4.19). Here, it is seen that 
small changes in kD; prove to have a large effect on _ the 
size of the small lobe at £ 99° and g° and 189° 7 EnUS 
yielding the largest error. A change of + 5% in kd effected 
the structure by up to 6 dB as seen from the behavior of the 
small lobe at + 90°. This change relates to a positioning 
error of about + 3 centimeters and presented a small problem 


with respect to exact positioning of the speakers. 


F. VALIDATION OF RESULTS 

The laboratory beam patterns prove to be in_= good 
agreement with theoretical beam patterns. Small uniform 
beam pattern amplitude differences are merely a consequence 
of the difference in normalization methods employed in the 
theoretical and experimental beam pattern processing. The 
theoretical program normalized the directivity to the sum of 
the amplitude weights, whereas the laboratory program 
normalized the beam pattern to the highest received voltage. 

An independent verification of the experimental program 
was performed by a thesis group under the auspices of 
Professor Don Walters of the Naval Postgraduate School 


Physics Department, who asked to use this experimental 
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system for related work. The laboratory system was placed 
in use by this research group in the selection of drivers 
for a 5 x 5 array currently under investigation and will be 


further utilized to establish beam patterns for this array. 
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TABLE [If 
ARRAY CONFIGURATIONS 
Run No/ Amplitude Spacing 


No. Preq Phase Tilt Shading from Center 
Elements : (Hz) +: (deg) : (deg) : (VAC) 1 (CM) 


ak 300 a 
l 38@ 45 
a 380 928 
i 9693 *) 
it 968 45 
at 968 98 
2 389 9,90 a 30 
2 300 8,180 ) 30 
2 962 8,0 @ 32 
2 9690 @,18@ ] 30 
Ps 962 8,2 98 34 
Z 962 @,18¢ 98 3a 
2 388 8,0 98 38 
Z 30a 8,189 98 32 
2 30a 8,20 45 30 
2 382 8,189@ 45 328 
2 962 0,9 45 328 
2 960 8,188 45 3a 
4 968 a,0,9,9 a 30 
4 30a 8,0,8,@ 8 38 
4 964 0,188,8,188 ’] 3a 
4 308 @,188,8,1828 a 38 
4 969 a,0,8,8 38 30 
4 382 0,9,3,2@ 308 30 
Quadll/4 300 8,188,3,188 30 3a 
Quad1l2/4 3a¢a 0,188,8,188 68 30 
Quadl4/4 96¢4 8,8,08,0 68 38 
Quad15/4 300 0,3,8,a 60 30 
Quadl6/4 960 0,9,0,9 90 30 
Quad17/4 96¢ 8,188,8,188 99 3a 
Quad18/4 968 0,188,8,20 a 30,30,38,9 
Triads 962 8,8,@ @ 3a 
Triz/s3 300 a,0,9 ) 30 
Tri 373 308 8,180,2 a 30 
Tri4d73 9624 @,180,2 ’) 30 
Tero 3 968 8,8,188 a 30 
Tri673 968 0,9,188 38 30 
Teva 3 9628 8,188,8 30 30 
TLeTe7 3 968 3,2,2 38 30 
Trvg973 96¢ 8,3,2a 30 30,15,3¢ 
Tri le73 964 9,3,3 a 3@.15,38 
Trilts 960 @,189,90 ’] 30,15,39 
Tries 3 962 6,189,9 30 30,15, 38 
TreiLol7 3 968 @,180,2 6@ 30,15,28 
Tril4/3 96@ 8,8,a 62 3¢@,15,3¢4 
Tritoss 3a 0,0,8 68 40°.15,30 
Tril6/73 968 0,0,2 98 30,15,30 
Trit77S 960 a,180,9 98 30,15,38 
Trpie7s 960 0,0,188 90 30,15,38 
Trv2o7 3 969 0,8,20 Q ee eres 30,9,38 
Tr12173 960 @,180,2 *] degial 30,3,390 
Tee227 3 960 a,08,2 ’] per a | 38,9,38 
Toiz2a7 3 968 8,180,2 ’] ae aes | 30,0,30 
Tete a7 968 0,8,189 ’] 1,272 30,0,30 
Tri257 3 968 0,23,2 f] pall 30,308,380 
Tri26/7 3 969 8,8,9 3a jy ek 38,30,38 
Trias 3 968 8,0,a 628 Lee a) 38,308,328 
Trizesa 962 8,8,8 92 5 Geer an | 38,308,380 
Trigd7a 968 8,8,@0 90 Lipa 38,38,38 
Tri3l7 3 968 8,188,28 92 pig We 38,308,380 
Trij2/3 968 8,8,9 a pe ee | 28,208,280 
Trio 962 8,188,29 a Letra 28,290,280 
Tevsa/ 4 960 9,180,28 a 2 20,108,298 
Toise/s 962 0,a,2a a iglespie 20,108,220 
Tris7 so 969 0,2,98 Q 5 eee eee 1aei8, 12 
Tri38/735 969 9,180,90 a Loe 18,108,120 
Rr 3o f5 969 a,189,2 92 5 re oI 20,380,280 
Tr44073 962 0,180,90 99 Deeg k “ley eos ad oo 
Tri4i73 964 0,08,2a 928 Lodged 1729517. oi 
Tri435/3 962 8,8,8 9a lo 2e 262526. 0 2050 
Tri44/3 962 8,188,28 98 1,2,1. 26.8,26.0,20.¢c 


Note: (1) Additional data runs were made in some 
of the same configurations with fiberglass 
wrapped around the centerpiece to identify 
the effects of scattering. 


(2) Additional data runs were made in some 


of the same configurations to verify 
repeatability. 
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Ve. CONCLUSIONS AND RECOMMENDATIONS 


The doublet forms the basis for countless sonar and 
radar arrays currently utilized in both military and 
Civilian applications. It is the ability to control the 
amplitudes and phases of the element of an array which 
allows the tailoring of the array radiation pattern to 
achieve the desired element of interest and promotes’ the 
maximization of array performance by increasing signal, 
decreasing noise or both. 

The ability to rapidly measure the beam pattern of an 
array with complete confidence in the results is the basis 
behind this’) study. A laboratory system was developed to 
allow collection of a 3698 point beam pattern ina 45 to 59 
minute period. Software is easily modified to allow for 
collection of more or less than 3698 data points. One of the 
major advantages of the computer-controlled method of 
obtaining a beam pattern is the relatively small amount of 
time required to achieve an accurate beam pattern based on a 
large data base. 

Recommendations for improving the current laboratory 
system are based on reducing electrical noise of the 
rotation motor and reducing the size and weight of the 
array. First, the noise of the rotation motor prohibits 


continuous rotation while reading the resistance on the 
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potentiometer. Since a CW signal does not require broadband 
detection, filtering could be accomplished by use of a lock- 
in, fast fourier transform or heterodyne filter. Second, the 
array weight should be reduced to lessen the drag on the 
rotation motor and to reduce the amount of coast-down time 
required for heavier arrays. Next, the array center joint 
should be reduced in size to lessen any shadowing effects 
from the large, metal center joint. (Scattering effects 
seemed to be small judging from the lack of any effect on 
the polar plots when this joint was wrapped in fiberglass.) 
Also, driver mounts need to be designed so that the center 
of the driver face is directly on the axis. The 2 
centimeter offset utilized for this investigation was 
minimized by making all readings in the far field. Last, 
greater distance between the driver and microphone would 
also help to lessen the 1 dB variation seen when an array 
was in an endfire orientation. This only affected the depth 


of the nulls. 
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APPENDIX A 


EQUIPMENT DESCRIPTION 


1. ANECHOIC CHAMBER 

All experimentation wasS conducted in the anechoic 
chamber, Room 819, Spanagel Hall, Naval Postgraduate School. 
The anechoic chamber is a 99% echo free enclosure for sounds 
above 1808 Hertz. Acoustics research can be conducted on 
sound sources, receivers and scatterers with minimal 
interference from Surface reflections or external noises. 
The size of the outer room in which the anechoic’ chamber 
floats is 37 feet x 24 feet x 21 feet high. Noise 
transmission from outside is minimized by the floating room 
wiekanm a room construction. The outer 12 inch concrete- 
walled room is separated from the inner room of concrete- 
block sides and floor by a 2 inch Dlanket of fiberglass or 
cork. Forty inch fiberglass wedges are attached to all 
surfaces of the inner room. These are designed to trap and 
absorb sounds’ which are incident upon them. A suspended 
floor is constructed of a grid of 225 wire cables, each 


stretched at a tension of 1598 to 289 pounds. 


2. HP-86B PERSONAL COMPUTER 
The HP-86B computer, operating in Hewlett-Packard BASIC 


language, iS equipped with a 64 Kilobyte (K-byte) memory 
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with two additional 128 K-byte memory modules. The system 
proved more than sufficient for both program development and 
as an instrument controller. Instrument control by means of 
the computer was accomplished with the HP-IB' (Hewlett- 
Packard Interface Bus) which is a standard connection on 
most HP equipment. This facilitated instrument control for 


the various system components. 


3. HP 7479A GRAPHICS PLOTTER 

The HP 747@A Graphics Plotter is a vector plotter which 
utilizes HP's micro-grip drive technology, moving both the 
paper and the pen. Multicolor plotting is available, but 
there are only two pen positions so pauses must be used to 
change pen colors. Only one line of program code is 
required to shift the output to and from the plotter to the 


HP=SOB sek i 


4, HP-4192A LOW FREQUENCY IMPEDANCE ANALYZER 

The HP-4192A LF Impedance Analyzer is a fully automated 
instrument designed to measure a large range of impedance 
related characteristics. It can be used to measure eleven 
impedance parameters: absolute value of impedance Z, 
absolute value of admittance Y, phase angleQ, resistance R, 
reactance X, conductance G, susceptance B, inductance L, 
capacitance C, dissipation factor D, and quality factor Q. 


The HP-4192A was used to measure the frequence' response, 
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impedance, resistance and reactance of the University Sound 
ID-38@ driver units as discussed in section 18 of this 


appendix. 


5. HP-3421A DATA ACQUISITION/CONTROL UNIT 

The HP-3421A Data Acquisition/Control Unit (DAC) has a 
Pe> digit voltmeter, a 10 kHz counter and an =4HP-IB 
controller interface. The DAC was utilized as a computer 
controlled switching unit with channels 1 and 2 configured 
aS actuators. This provided a voltage path to control 
relays in the motor power and direction controller for 


rotation of the ID-398 speaker array, Fig. (A.l). 


6. HP-3478A MULTIMETER 

The HP-3478A Multimeter provides both resistance and 
voltage measurement capabilities with 3.5 to 5.5 digit 
resolution with + @.3 mV accuracy on the 100 Hz to 20 kHz 
scale. The multimeter was remotely controlled by the 
computer via the HP-IB. Voltage readings from the GR-1962 
microphone and resistance readings from the variable 
potentiometer on the array rotator being read by the 
multimeter and then sent to the HP-86B computer via the HP- 
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73 HP-3314A FUNCTION GENERATOR 
The HP-3314A Function Generator is a multimode function 


generator which provides” sine, Square, and triangular 
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Data Acquisition / Control Unit 
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to HP Multimeter 
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Power Suppl 
ee to Motor Power & 


Direction Controller 


Figure A.1. Data acquisition connection to motor-power 
controller. 
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functions as well as any desired arbitrary wave form from 
mecouencies Of 8.09IeHze= to 19.999 MHz. It was operated in 
the free run mode providing a continuous 1 Vac sine wave to 


Hp-467A Amplifier's inputs. 


8. HP-465A AMPLIFIER 

The HP-465A Amplifier is a general purpose preamplifier 
and impedance converter from 18 megohms to 58 ohms. It has 
a selectable gain of 28 or 48 decibels and is stable over a 
frequency range of 5 Hz to 1 MHz. The amplifier was used to 


amplify the received signal from the GR-1962 microphone. 


9. HP-467A POWER AMPLIFIER 


The HP-467A Power Amplifier is a 18 watt peak power 
amplifier and a -20 to +20 vold dc power supply. It has a 


wide band width with low de drift from 5 Hz to 1 MHz. A 
variable gain control allows adjustments of gain between 
zero and 18 with resolution better than 9.1% of full power. 
Four HP-467A's, one for each ID-3@ driver unit, were 
utilized in adjusting power inputs from the function 


generator. 


18. UNIVERSITY SOUND ID-38 DRIVER UNITS 

The University Sound ID-3@'s are medium power horn 
driver units designed specifically for wide midrange 
response. Their tropicalized voice coils are 2 inches in 


diameter with linen base molded phenolic diaphragms with a 


zea 


minimal 16 ohm impedance. The frequency response is 85 Hz 


to 7.5 Khz. They were used as the sound sources. 


ll. GENERAL RADIO TYPE 1962 MICROPHONE 

The GR-~-1962 microphone is a 9.5 inch diameter electret 
condenser, microphone, providing a flat response at zero 
decibels between 280 Hz and 15 Khz. It experiences'7 rapid 
drop off above 15 Khz and is nonresponsive above 21 Khz. A 
General Radio Model 1568~-P62 inline power’ supply was 
connected to the microphone 1568-P42% preamplifier to 


amplify the received signal. 


12. HP-6215A POWER SUPPLY 
The HP-6215A power supply provides zero to 38 volts dc 
or ®@ to 508 Mamps. It was used to provide switching voltage 


for the Motor power relay in the direction controller. 


13. MOTOR POWER AND DIRECTION CONTROLLER 

The controller for the array rotator controlled power to 
the rotator and also the direction of rotation. With leads 
to the rotator assembly potentiometer, it also provided 
resistances which corresponded to a given angle (8 to 19 
Kohms translated into 8 to 368 degrees). A +5 vVOltUgiGe 
signal to IC-l via channel 1 of the DAC provided on/off 
power to the rotator and +5 Volts to IC~-2 via channel zero 
of the DAC provided clockwise rotation of the driver 


assembly unit (Figure A.2). [Ref. 8] 


eZ 











Dip Relay Input and Output 


AC high 
AC low 
+3 V from DAC Chi 
Common 

+5 V from DAC Ch O 
Ohm sensing High 
Ohm sensing Low 


ee ae 


© Cris 

© 4 13 

6 11 

Q 8 ; 
b) 360° / 10 kM Potentiometer c) 16 Pin DIP Relay 


Figure A.2. Circuit diagram for modified motor-power / direction 
controller (Ref. 9., p. 57]. 
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APPENDIX C 


Experimental and Theoretical Plots 
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